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ABSTRACT 

If dark halos are composed of dense gas clouds, as has recently been inferred, then 
collisions between clouds lead to galaxy evolution. Collisions introduce a core in an 
initially singular dark matter distribution, and can thus help to reconcile scale-free 
initial conditions - such as are found in simulations - with observed halos, which have 
cores. A pseudo-Tully-Fisher relation, between halo circular speed and visible mass 
(not luminosity), emerges naturally from the model: M V i 8 oc y 7 / 2 . 

Published data conform astonishingly well to this theoretical prediction. For our 
sample of galaxies, the mass- velocity relationship has much less scatter than the Tully- 
Fisher relation, and holds as well for dwarf galaxies (where diffuse gas makes a sizeable 
contribution to the total visible mass) as it does for giants. It seems very likely that this 
visible-mass/velocity relationship is the underlying physical basis for the Tully-Fisher 
relation, and this discovery in turn suggests that the dark matter is both baryonic and 
collisional. 
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INTRODUCTION 

A great variety of dark matter candidates exist, motivated 
by diverse pieces of evidence, typically indirect (see, for 
example, the review by Trimble 1987). One such piece 
of evidence comes from the "Extreme Scattering Events" 
(ESEs: Fiedler et al. 1987); these are radio-wave lensing 
events caused by dense blobs of plasma crossing the line-of- 
sight. Walker & Wardle (1998a) presented a model in which 
these events are caused by ionised material associated with 
planetary-mass, molecular gas clouds in the Galactic halo. 
This model is good at explaining the ESE phenomenon, but 
carries with it the implication that most of the mass of the 
Galaxy is in this cold, dense form. If the Galactic dark mat- 
ter is really in cold gas clouds - as, in fact, has been proposed 
previously by a number of authors (Pfenniger, Combes & 
Martinet 1994; de Paolis et al 1995; Gerhard & Silk 1996) - 
then consistency with a variety of data requires that these 
clouds satisfy several constraints (Gerhard & Silk 1996). 
Foremost amongst these is the requirement that collisions 
between clouds should not entirely deplete the halo of its 
dark content. In this paper we use the simplest halo model, 
an isothermal sphere, to show how cloud-cloud collisions lead 
directly to structural changes in the halo (§2), and how this 
process leads to predictable visible galaxy masses (§3); how- 
ever, no attempt is made to predict the form which the vis- 
ible mass should assume (stars vs. diffuse gas). 

2 ISOTHERMAL HALOS 

In a halo with one-dimensional velocity dispersion a the typ- 
ical relative speed of a pair of clouds is y/Ea. Essentially all 



collisions are highly supersonic and even glancing impacts - 
i.e. those with impact parameter roughly equal to twice the 
cloud radius - are expected to unbind the clouds, with the 
result that the collision products become visible as diffuse 
gas (which may subsequently be transformed into stars). In 
this circumstance the cross-section for disruptive collisions 
between clouds is just four times the geometric cross-section 
of a single cloud. Recognising that each collision disrupts two 
clouds, we see that the collision rate can be written as 



dlog c P 
dt 



id 



where p is the halo density, and E is the mean surface-density 
of the individual clouds. A more rigourous treatment, in 
which one integrates the collision rate over the velocity dis- 
tribution of the clouds, yields a numerical factor (32/^/n) 
which differs by less than 10% from that in eq. 1. 

Collisions occur preferentially between pairs of clouds 
having anti-parallel velocities, so that the collision products 
have a much smaller velocity dispersion than the parent 
clouds (by a factor of v8), and subsequently undergo in- 
fall in the gravitational potential. This infall modifies the 
potential, which in turn leads to evolution of the dark halo 
density beyond that described by equation 1. (This is "halo 
compression" — see Blumenthal et al 1986.) We shall not 
attempt a self-consistent treatment but, rather, we neglect 
the evolution of the gravitational potential. This allows us to 
estimate the dark halo density evolution straightforwardly 
by integrating equation 1; the result is 



p(r,t) = 



l + t/t c (r)' 



(2) 
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Figure 1. Data for a sample of 15 galaxies with well determined rotation curves (see Table 1), showing the total visible mass as a function 
of galaxy circular speed (i.e. M v i a [V]). Contributions to M v i s come from stars - both bulge and disk (determined from "maximum disk" 
decompositions) - and diffuse gas. The plotted symbol size is arbitrarily chosen. The line shows the prediction of equation 9, with cloud 
surface density (E) treated as a free parameter. 



where the characteristic time t c (r) = Y,/8\/Hap(r, 0). If we 
now model the initial distribution of dark matter as a sin- 
gular isothermal sphere, i.e. p(r, 0) = a 1 /2tyGt 2 , for radius 
r, then equation (2) implies 



p(r,t) = 



2nG{r 2 +ri)' 



where 



2 

r r = 



AV6 a 3 t 
n GE' 



(3) 



(4) 



That is, at t > the halo is a non-singular isothermal sphere 
with core radius r c . We can recast equation 4 in terms of 
the native circular speed, V, for the halo, using V = \/2a. 
Taking the proposed mean column for individual ESE clouds 
(Walker & Wardle 1998a; see also §4) of TV ~ 10 25 cm" 2 , 
i.e. E ~ 30 g cm -2 , together with t = 10 Gyr, we obtain 



4.2V-/00 5 



kpc, 



(5) 



where V = 100 V100 kms -1 . In consequence, at the present 
epoch we should expect galaxy halos to possess modest 
cores, if those halos virialised at redshifts 2 > 1 (i.e. 
~ 10 Gyr ago). These core radii ought, in principle, to be 
measurable from rotation curves of spirals. 

In practice the task of measuring a core in the den- 
sity profile of a dark halo is made difficult by the visible 
galaxy, which itself contributes to the rotation curve, and 
one is obliged to perform a disk-halo decomposition — a 
process which engenders a number of uncertainties. Further- 



more, the visible galaxy can alter the gravitational potential 
sufficiently that the halo density distribution changes in re- 
sponse: "halo compression"; see Blumenthal et al. (1986). 
Notice that this effect is not accounted for in equation 5. 
These difficulties are mitigated to a large extent if the visi- 
ble galaxy is so puny that it acts only as a kinematic tracer, 
allowing the rotation curve to be measured but not playing 
any significant role in determining its form. Such galaxies ex- 
ist and have already been used to demonstrate that, at least 
in these cases, halos do possess cores (Moore 1994; Flores & 
Primack 1994). This is an important result because it ap- 
pears to conflict with simulations of halo formation using 
dissipationless dark matter: these exhibit singular density 
profiles (Dubinski & Carlberg 1991; Navarro, Frenk & White 
1996; for an alternative view see Kravtsov et al. 1998). 

Kormendy (1992) has derived scaling laws for real (late- 
type spiral) galaxy halos; his relationship between core ra- 
dius and halo circular speed is (for Ho — 75 kms" Mpc~ ) 



r c = 3.6 V100 kpc; 



(6) 



encouragingly close to equation 5. This scaling is less secure 
than the simple fact that cores exist in galaxy halos, and one 
should beware of undue emphasis on the numerical agree- 
ment between equations 5 and 6, as the mean cloud column 
density (E) is really only known to within a factor of a few 
(see also the comments in §4). Indeed in the following sec- 
tion we shall treat E as a free parameter, and the resulting 
best fit indicates a larger value of E than used here. 
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Figure 2. As figure 1 but for the stellar mass only, i.e. excluding the contribution of diffuse gas to the total visible mass; the theoretical 
line is identical to that in figure 1. It is clear that the dwarf galaxies, which have a high mass-fraction in diffuse gas, systematically depart 
from the relationship defined by the giants, if only stellar mass is included in the accounting. Compare this to figure 1, where a single 
relation is valid for the whole sample. 



Table 1. Galaxies with well determined rotation curves. Circular speeds are taken from table A.l (column 9), p94 of Brocils (1992); 
while masses are obtained from "maximum disk" rotation-curve decompositions — the sum of columns 7 (diffuse gas), 9 (disk stars) 
and 11 (bulge stars) of table 2, p244, in Brocils (1992). Only 15 galaxies are common to both tables. 
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(kms- 1 ) 


(10 10 Mq) 


NGC55 


86 


0.42 


NGC247 


108 


1.05 


NGC300 


90 


0.58 


UGC2259 


90 


0.40 


NGC1560 


74 


0.22 


NGC2403 


134 


1.82 


NGC2841 


294 


35.0 


NGC2903 


180 


4.82 


NGC3109 


64 


0.10 


NGC3198 


149 


3.80 


DD0154 


45 


0.05 


NGC5033 


220 


10.8 


DDO170 


60 


0.09 


NGC6503 


115 


1.05 


NGC7331 


238 


18.0 



3 A PHYSICAL BASIS FOR THE 
TULLY-FISHER RELATION 

The Tully-Fisher (TF) relation is an empirical result which 
connects the width of 21 cm line emission, AV, from a spiral 
galaxy with the galaxy's expected luminosity, L: 

LcxAV Q , (7) 

where a ~ 4 (Strauss & Willick 1995); the scatter in this 
result is quite small. Unfortunately there is no sound the- 
oretical basis for the TF relation. Attempts to derive it on 
the basis of the virial theorem, plus the assumption of fixed 
surface- brightness for all galaxies ("Freeman's law"), fall 
foul of the fact that TF holds also for very low surface- 
brightness galaxies (Zwaan et al 1995; McGaugh & de Blok 
1998). Put another way, we expect that AV ~ 2V, so that 
the kinematics reflect the properties of the dark halo, but L 
is a manifestation of the stellar component of the galaxy, and 
no direct coupling between these two is expected a priori. In- 



deed, since total mass-to-light ratios can vary by more than 
an order of magnitude amongst the population of observed 
galaxies, one tends to think of the visible and dark compo- 
nents as almost independent. However, the collisional pro- 
cess we have described in §2 converts dark matter into visible 
forms, and thus creates a close coupling between the dark 
halo and the visible galaxy. Such a coupling holds promise 
for explaining the Tully-Fisher relation; this appealing at- 
tribute of baryonic models has been recognised previously 
(Pfenniger, Combes & Martinet 1994; Gerhard & Silk 1996). 

We can compute the total visible mass, M v u, in the 
form of stars and gas (without actually saying anything 
about their relative proportions), from 

M vi3 (t)= / drAnr 2 [p(r,0)- p{r,t)]. (8) 
Jo 

With the evolving isothermal halo described by equation 3, 
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this gives 

2 

M vis = ^~r c . (9) 

Notice that M v i s is perfectly well defined, despite the total 
mass of the dark halo diverging at large radii. If we adopt 
the reasonable supposition of a roughly similar stellar mass- 
to-light ratio for all spirals then, provided that diffuse gas 
makes a negligible contribution to M v i s , we expect L oc V 3 ' 5 , 
which is close to equation 7. In the usual form, however, TF 
is a relation between global 21 cm line- width, AV, and lumi- 
nosity. Moreover, various exponents a ~ 4 are observed, so 
it is not immediately obvious that our theory is in agreement 
with the data. 

Spectral line imaging (notably 21 cm imaging) gives 
detailed information on the velocity field of a galaxy, allow- 
ing accurate determination of the rotation curve (hence V), 
rather than just the global property AV which can be de- 
termined with a single dish radio telescope. Broeils (1992, 
chapter 4, appendix A) has investigated the L[AV] and L[V] 
relationships for a sample of 21 galaxies with well deter- 
mined rotation curves. He finds that blue luminosity is more 
tightly correlated with V than with AV; the former displays 
a scatter of 0.22, and the latter 0.28, in log 10 L. (More pre- 
cisely, Broeils used the circular speed in the flat part of the 
rotation curve, which we take to be a good estimator for V .) 
Furthermore, the relation he derives is L oc V 3 ' 4 , which is 
very similar to the form we predict. 

We could go further and compare the normalisation of 
this result with that of our own theory, but in doing so 
we would be forced to introduce an additional quantity, the 
value of the stellar mass-to-light ratio, which is a priori un- 
known. This can be avoided if we compare M V i S directly 
with mass estimates derived from rotation curve decompo- 
sitions which, in effect, measure the stellar mass-to-light ra- 
tio for each galaxy. Furthermore, by utilising mass, rather 
than luminosity, we expect better agreement with our the- 
ory because we can include the diffuse gas content. Note that 
diffuse gas is often a substantial fraction of the total visible 
mass of dwarf galaxies; in the (extreme) case of DD0154 it 
amounts to 80% of the visible mass (Carignan & Freeman 
1988). While Broeils (1992) did not investigate the possi- 
bility of any correlation between M v i s and V, he did give 
maximum disk decompositions of the rotation-curves for 15 
of the 21 galaxies studied; his results are summarised in 
table 1, and plotted in figure 1. Also shown in figure 1 is 
the relation given by our equation 9, with the cloud surface 
density (E) treated as a free parameter; the fit is manifestly 
good. 

The very first point which must be made about this re- 
sult is that it is not simply an artefact of (maximum-disk) 
rotation curve decomposition. This point is clear because 
rotation speed measures M/R, not M, so even if the visible 
mass were the dominant component we would still be obliged 
to explain the existence of a characteristic length scale hav- 
ing a form similar to that of equation 4. However, in the 
flat part of the rotation curve - which has been used to es- 
timate V - the visible matter contributes less than the dark 
matter in every case in table 1. On average the dark matter 
contributes three quarters of the total, for this sample, and 
for some galaxies the visible contribution is negligible. The 
hypothesis that figure 1 simply reflects the rotation curve 



decomposition procedure is therefore untenable. 
The explicit numerical form of our fit is 

M vls = 7.0 x 10 9 V^ 2 M , (10) 

corresponding to E ~ 140 gem -2 (for t = 10 Gyr), and 

r c = 1.9V 1 3 / 2 kpc. (11) 

We will not attempt to give a figure of merit for the fit 
quality, as the uncertainties associated with the data in- 
volve systematic uncertainties in rotation-curve decomposi- 
tion, and these are hard to quantify. We can, however, mea- 
sure the scatter of the data about the theoretical prediction: 
the root-mean-square deviation is 0.084 in log 10 M v i s . By 
comparison the scatter in the Tully-Fisher relation for this 
sample is almost a factor of 3 larger (0.24 in log 10 Ls - 
slightly bigger than the scatter of 0.22 for the full sample 
of 21 galaxies studied by Broeils). Notice that only a part 
of this difference can be accounted for by uncertainties in 
galaxy distances, as the stellar mass inferred from rotation 
curve decomposition is proportional to distance, while lumi- 
nosity and the measured gas mass are both proportional to 
(distance) 2 . We also note that if the mass contributed by dif- 
fuse gas is neglected then the dwarf galaxies systematically 
depart from the relation defined by the giants; this point is 
graphically illustrated in figure 2, where stellar mass is plot- 
ted as a function of halo circular speed. These facts oblige us 
to conclude that the fundamental connection is between total 
visible mass and halo circular speed, with the Tully-Fisher 
relation emerging as an approximation which is valid when 
most of the visible mass is in stellar form. 

To add a little more weight to this conclusion, we em- 
phasise that our sample of galaxies is very heterogeneous: 
it spans the entire size spectrum from dwarfs to giants; it 
includes low surface-brightness objects; most importantly, 
perhaps, it includes cases in which the visible galaxy makes 
a negligible contribution to the rotation curve. This last 
point is crucial as it requires that the observed M„j s [V] re- 
lation be interpreted as M vis being determined by V, and 
not the other way around. That is, purely from an observa- 
tional perspective we can assert that the visible mass content 
of a galaxy is determined by the velocity dispersion of the 
dark matter halo. Our theory shows why this ought to be 
so, and it follows that these data support the model of a 
baryonic dark halo, with clouds of typical surface density 
E ~ 140 gem" 2 . 

4 DISCUSSION 

The evolution implicit in equation 9 could, in principle, be 
used to further test the theory we have presented, but in 
practice this will be difficult. In particular the requisite sen- 
sitivity and angular resolution, for determining 21 cm rota- 
tion curves of normal galaxies at z ~ 1, are both beyond 
the reach of current instrumentation. Some studies have al- 
ready been made of Tully-Fisher-type correlations of galax- 
ies at z ~ 1, based on optical data alone (e.g. Vogt et al. 
1997). These show mild evolution in a sense opposite to that 
expected in our theory (if a constant mass-to-light ratio is 
assumed), and we must suppose that this is a result of stel- 
lar populations being younger at earlier epochs. Notice that 
there are serious consequences for galaxy distance estimates 
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based on a local Tully-Fisher relation if the actual TF rela- 
tion evolves with look-back time. 

It is well known that the infrared TF relations show less 
scatter than their visible counterparts (Aaronson, Huchra & 
Mould 1979), because infrared photometry is less sensitive 
to the young, hot stars which tend to dominate the luminos- 
ity (but not the total stellar mass) of the galaxy. Extinction 
corrections are also smaller in the infrared. This argues that 
rotation-curve decompositions based on infrared photome- 
try ought to be better than those based on visible photom- 
etry. In addition, there is no particular reason to suppose 
that "maximum disk" rotation curve decompositions give 
the best estimates of the stellar mass-to-light ratio. Thus 
with a study that is specifically aimed at testing equation 
9, it may be possible to reduce the scatter seen in figure 1 
below its already small value. Ideally one would like to work 
with a sample of galaxies from a single cluster, thereby re- 
ducing the dispersion contributed by distance uncertainties; 
distance errors probably dominate the currently observed 
scatter of 0.084 dex in the M„j s [V] relation. 

Any cluster of galaxies which is well approximated by 
an isothermal distribution of dark matter should fit into 
the scheme we have outlined, simply by an appropriate 
choice of velocity dispersion. (We estimate from the evapo- 
ration constraints given by Gerhard & Silk [1996] that the 
clouds can survive for roughly a Hubble time in the en- 
vironment of a rich cluster of galaxies.) For example, tak- 
ing a = 10 a kms -1 equation 10 predicts a visible mass of 
M V is — 7 x 10 13 Mq , broadly consistent with the data (Jones 
& Forman 1984), and equation 11 gives a dark matter core 
radius of r c = 100 kpc. This core radius is rather larger than 
indicated by analysis of cluster lensing data (Miralda-Escude 
1993; Flores & Primack 1994) but we must bear in mind 
that lensing measures the total surface density, not just the 
dark component. For galaxies the greatest visible contribu- 
tion typically comes from stars, whose presence is relatively 
straightforward to quantify, but for rich clusters it is usually 
the hot intra-cluster medium which dominates, and here the 
inferred mass distribution is much more model dependent. 
(Notice, again, that the theory we have presented is only par- 
tially predictive in that it gives the total visible mass, but 
does not tell us whether this ought to be in stars or diffuse 
gas.) In particular we note (i) it is thought that in "cooling 
flow" clusters a large amount of gas accumulates in some 
(unknown) form at the centre of the cluster (Fabian 1994), 
and (ii) there are tentative detections of huge amounts of 
warm (EUV emitting) gas in some clusters (Mittaz, Lieu & 
Lockman 1998). 

In the foregoing discussion we have given no consider- 
ation to the material properties of the cloud collision prod- 
ucts, being content with the notion that this stuff becomes 
part of the visible pool. A basic analysis of the effects of 
the shock (Walker & Wardle 1998b) indicates that, for the 
Galactic halo, the result of a collision will typically be atomic 
gas, even though the clouds are initially molecular. This sug- 
gests a possible connection with the "High Velocity Clouds" 
(Wakker & van Woerden 1997) which are seen (almost ex- 
clusively) in 21 cm emission: some HVCs may simply be 
material from recent dark cloud collisions. This gas is ex- 
pected to have such a low column density that it will be 
stopped by the diffuse interstellar medium, at the first en- 
counter with the Galactic plane, thereby contributing to the 



assembly of the gaseous disk and, subsequently, the Milky 
Way. 

As well as being relevant to the Tully-Fisher relation, 
the model we have presented is germane to the "Disk-Halo 
Conspiracy" (Sancisi & van Albada 1987). The conspiracy is 
so-called because it seems puzzling that (giant) spiral galaxy 
rotation curves should be as flat as they are observed to be, 
given that the acceleration is predominantly due to stars 
at small galactocentric radii, and due to dark matter at 
large radii. In the model we have presented, however, the 
conspiracy is no surprise; rather it is an innate feature, as 
dark matter is converted to visible forms by cloud-cloud col- 
lisions. With the model of §2 we cannot sensibly compute 
rotation curves - because we have neglected evolution of the 
gravitational potential - so these statements are necessarily 
qualitative. A useful quantitative treatment would require 
a self-consistent description of the evolution; more careful 
consideration of the initial conditions (dark matter density 
profile) would also be appropriate. 

The principal difficulty with our model is that the fit 
shown in figure 1 requires E ~ 140 gem -2 (for t = 10 Gyr), 
whereas Kormendy (1992) measures halo core radii which 
suggest a smaller cloud surface density, E ~ 40 g cm -2 . Both 
of these values are higher than the original estimate (Walker 
& Wardle 1998a) based on scintillation data, and incompat- 
ible with each other. Kormendy (1992) cautions that his 
derived halo scaling relationships are quite uncertain, and 
it is evident from his figure 3 that the scatter in the mea- 
surements is much larger than in our figure 1, so the value 
S = 140 gem -2 is to be preferred. One possible resolution is 
that our theory is simply too crude. For example, in the case 
of spiral galaxies the fact that the visible matter has angular 
momentum requires that the dark matter should also have 
some, yet we are modelling halos as non-rotating. Perhaps 
a more refined theory would yield consistent estimates of 
halo core radius and visible mass? We have made a prelimi- 
nary investigation of this possibility, using Toomre's (1982) 
analytic models of rotating isothermal halos. These models 
span the entire range from isothermal spheres to cold, ro- 
tationally supported disks. We find that that the prediction 
M vi3 oc V 1 1" 1 holds for this whole set of models, with a co- 
efficient of proportionality which is only weakly dependent 
on the degree of rotational support that the halo possesses. 
Similarly the inferred core radius is only weakly dependent 
on rotation, so it seems unlikely that the discrepancy can 
be explained in this way. The fact that the predicted M v i s 
is only a weak function of the degree of halo rotation, for a 
fixed value of V does, however, lend a certain robustness to 
the model we have presented. 

Finally we note that the preferred cloud column density 
deduced in this paper, i.e. E ~ 140 gem -2 , is substantially 
larger than the value E ~ 30 g cm -2 anticipated on the basis 
of the ESE data (Walker & Wardle 1998a). The latter value 
was derived from the sky covering fraction, r, of the ESE 
clouds, together with the requirement that the cold clouds 
not contribute more than the dynamically determined mass 
for the Galactic dark halo. Using the density distribution of 
equation 3, we find that at high Galactic latitude, 



at the Solar circle, where E = lOOEioo gem 2 . Evidently if 



6 Walker 



Eioo = 1-4, then the dynamical constraints limit the cover- 
ing fraction of the clouds to be r ^ 2 x 10 -4 . The observed 
ESE cloud covering fraction is highly uncertain because 
what actually constitutes an ESE is only very loosely de- 
fined. Selecting all periods of "unusual variability", Fiedler 
et al (1994) reported T t, s ~ 5 x 10 -3 . On the other hand, 
if we consider events whose light-curves are reminiscent 
of the calculations of Walker & Wardle (1998a), then we 
would include only the events for the sources 0954+658 and 
1749+096, leading to an estimate Tobs ~ 5 x 10~ 4 ; this is 
still 2.5 times larger than we expect for E = 140 gem -2 . 

To understand how this discrepancy might arise we need 
to recognise that the size of an ESE cloud relates to the inner 
boundary of the ionised wind, as it is the free electrons which 
comprise the radio- wave "lens". By contrast the measure of 
E we have given in this paper is a characteristic surface 
density for the underlying hydrostatic cloud; the wind is 
irrelevant here as it contains a negligible fraction of the mass 
of the cloud. Walker & Wardle (1998a) explicitly assumed 
that the inner boundary of the ionised wind coincided with 
the surface of the hydrostatic cloud. However, if the inner 
boundary of the ionised wind is actually the outer boundary 
of a neutral wind, then the ESE covering fraction will be 
substantially larger than the estimate given by equation 12, 
as is observed to be the case. It is plausible that a neutral 
wind should underly the ionised wind, but this possibility 
requires careful consideration in its own right and will not 
be pursued here. Further constraints on the column density 
of the cold clouds are considered by Draine (1998), based on 
optical refraction by the neutral gas ("gas lensing"). These 
constraints are compatible with E ~ 140 gem -2 provided 
that the internal density profile of the clouds is not strongly 
centrally concentrated; convective polytropes, for example, 
are acceptable. 



5 CONCLUSIONS 

Modelling galaxy halos as isothermal spheres composed of 
collisional, baryonic dark matter leads us to expect non- 
singular dark halo density distributions with a predictable 
visible mass content. Both of these expectations are borne 
out in the data, with the latter result very likely being the 
fundamental basis for the Tully-Fisher relation. Although 
simple, the theory is remarkably good at predicting (spiral) 
galaxy masses across the whole size spectrum from dwarfs to 
giants. The evident success of this description of galaxy evo- 
lution gives strong support to the notion that galaxy halos 
are composed of a vast number of cold, dense, planetary- 
mass gas clouds. 
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